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Atomistic studies are carried out to investigate the formation and evolution of defects during 
nanoindentation of a gold crystal. The results in this theoretical study complement the experimental 
investigations [J. D. Kiely and J. E. Houston, Phys. Rev. B 57, 12588 (1998)] extremely well. The 
defects are produced by a three step mechanism involving nucleation, glide and reaction of Shockley 
partials on the {111} slip planes noncoplanar with the indented surface. We have observed that slip is 
in the directions along which the resolved shear stress has reached the critical value of approximately 
2 GPa. The first yield occurs when the shear stresses reach this critical value on all the {111} planes 
involved in the formation of the defect. The phenomenon of strain hardening is observed due to the 
sessile stair-rods produced by the zipping of the partials. The dislocation locks produced during the 
second yield give rise to permanent deformation after retraction. 

PACS numbers: 62.20.Fe, 62.20.Qp 



I. INTRODUCTION 

Understanding the detailed mechanics of material deformation is a fundamental challenge in materials science. Ju. 
metals, the defect structures produced during deformation influence the material properties and behavior critically.Em 
The formation and evolution of such structures have their basis in atomistic processes and the study of these nanoscale 
phenomena is paramount to the understanding of macroscopic phenomena such as fracture, friction, strain hardening 
and adhesion. The results of such research will also greatly facilitate the design of novel materials with desired 
properties. These insights into material behavior can be exploited to create desirecL-dislocation patterns which can 
then be etched in a controlled manner to fabricate nanopatterns and nanostructurcs.u 

Nanoindentation experiments, with the advent of scanning probe microscopes and advances m-.iiidentation tech- 
niques, are capable of experimentally probing material properties and phenomena at the nanoscale .Qtffi At these atomic 
length scales, the continuum models of deformation do not perform well and atomistic methods need to be considered 
to investigate the nanoscale deformation behavior. Advances in computational capability and high performance tech- 
niques have enabled researchers .^investigate nanoindentation studies of comparable length scales theoretically using 
molecular dynamics simulations.l3u The experiments, for most part, have emphasized quantitative investigation of me- 
chanical properties by measuring the force displacement curves, and the theoretical computer simulations have been 
targeted at studying the atomistic processes involved in plastic deformation during indentation experiments. The pri- 
mary goal of such studies is to complement the experimental findings with theoretical investigations in understanding 
the mechanisms of plastic deformation in materials. 

In this paper, we present results of atomistic studies of nanoindentation of a passivated gold surface. The objective 
of this work is to study the atomistic processes responsible for plastic yield duting the initial stages of indentation 
and explain the experimentally observed yield phenomena and defect structures .0 

II. METHODOLOGY 

The objective of this atomistic study is to investigate the defect nucleation during nanoindentation of a passivated 
Au(OOl) surface and study the mechanisms leading to plastic deformation. The atomic configuration of the system 
studied is illustrated in Fig. The gold substrate is modeled as a slab (122 A x 122 A x 50 A ) containing 46400 
atoms with periodic boundary conditions parallel to the surface. The orientation of the slab is such that the directions 
[100], [010] and [001] coincide with x, y and z. The bottom layer is fully constrained and the substrate size is sufficiently 
large to eliminate the finite size effects. The indenter is an assemblage of atoms in diamondoid cubic lattice arranged 
as a truncated pyramid with exposed (111) facets and a 15 A x 15 A (001) square indenting face. The indenter is 
oriented such that the edges of the indenting face arc in [110] and [110] directions with respect to the gold crystal. 

We have employed the quantum Sutton-Chen (Q-SC)El potential to model the gold atoms. This formulation includes 
the quantam corrections to better predict mechanical properties, and retains the simplicity of the original Sutton-Chen 
potentiaM to facilitate the understanding of the underlying physics of various processes. The indenter is modeled as 
a rigid body and the indenter-surface interactions are purely repulsive, V(r) = e(r/p)~ 12 with e = 25 meV and p = 
3 A , to eliminate the adhesive interactions and mimic the passivation of the gold surface in experiments. 



FIG. 1: Atomic configuration of the indenter and the gold 
substrate. • - indenter; O - dynamic gold substrate; O - 
temperature control region; • - fully constrained boundary. 
The indenting face is a square with edges along the [110] and 
[1 10] directions of the gold crystal. 



We have used an extended version of the parallel MD package DL_POLYM to perform the calculations. The 
dynamics of the substrate is evaluated by integrating the Newtonian equations of motion using Verlet-leapfrog method 
with a timestep of 1 fs. The gold substrate is equilibrated to its minimum energy configuration at 300 K and the 
indentation is simulated by advancing the indenter atoms by 0.0005 A at every timestep, giving the indenter a 
velocity of 50 m/s. The temperature is regulated by periodically scaling the velocities of the atoms of the deepest 
non-constrained region of the substrate, away from the contact region to minimize the interference of the temperature 
control mechanism with the normal energy flow processes that occur in the contact region. ■— ■ 

To understand the mechanics of plastic deformation during indentation, atomic stress tensorlia er is used to study 
the distribution of stresses. The von Mises shear stress, -y/jj, proportional to the square root of the distortion energy, 
is an indicator of the onset of plastic yieldingcJ as proposed by von Mises. The von Mises shear stress is given by the 
square root of the second invariant of the deviatoric stress, J2, which is defined as 



where Tr denotes the trace of a matrix, I is the unit matrix and p is the local hydrostatic pressure. 

In metals, plastic deformation occurs by the glide of dislocations on the slip planes. In ordermto identify and 
characterize the dislocations being nucleated during indentation we employ the slip vector analysis,El which provides 
information on the Burgers vectors of dislocations. The slip vector is defined as, 



where, n s is the number of slipped neighbors /?, of the reference atom, and rf and r p B are the vector differences of 
atom (3 and the reference atom positions at times t and 8, respectively. The slip vector given by the above expression 
represents the Burgers vector of slip between the plane containing atom a and its adjacent atomic planes, in the time 
interval [0,t\. However, this is true only in the case of single slip, where the reference atom is contained by only one 
slip plane. In the event of multiple slip, where the atom is contained by two planes undergoing slip simultaneously, the 
Burgers vector is different from the slip vector given above. In any event, the slip vector will have a large magnitude 
for inhomogeneous, non-affine deformation near the atom and can be used to identify slipped regions. 

The strains, induced by indentation are studied by evaluating the atomic strain tensor as formulated by Horstemeyer 
and Baskcso This formulation is based on the deformation gradient for a material employing many-bodied potential. 
The atomic Lagrangian Green strain tensor E, used in this study is given by, 



J 2 =±Tr [(*- pi) ■ (ct - plf] , 
p=-iTr(<r). 
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where, F is the deformation gradient, m is the number of nearest neighbors /3 of the reference atom, ig,rf have 
the same meaning as above and <g> represents tensorial product. E will then quantify the strain experienced by the 
reference atom in the time interval [9,t\. 

To investigate the mechanisms of dislocation nucleation and glide on the slip planes, we study the resolved shear 
stresses on the slip planes along the Burgers vectors given by the slip vector analysis. The resolved shear stress r on 
a plane with normal n along the direction of slip b is given by, 



T ( n )[b] =b cr n (5) 

Schmid lawEl states that a slip system is activated when the resolved shear stress on that system reaches a critical 
value called the critical resolved shear stress (CRSS). 



III. RESULTS AND DISCUSSION 



A. Indentation 

The force versus displacement curve for the initial stages of indentation is shown in Fig. ||. The force F z is 
calculated as the sum total of the forces exerted on the indentcr atoms by the substrate and the displacement z sep 
is the separation between the indenter apex and the surface of the substrate before indentation. Initially, the force 
curve displays elastic behavior until the force decreases abruptly at the first yield point, marked as (1) in Fig. [| This 
phenomenon is associated with the nucleation of a plastic event to parlijlly_xclieve the elastic stress in tbe_cpntact 
region. This observation is in excellent agreement with other theoreticaflaHEZEa and experimental results. D'ErEI Upon 
further indentation the force begins to rise again, displaying yet another region of elastic behavior, until the substrate 
undergoes a second yield event (2) in Fig. g. It is interesting to see that the force curve has a higher slope in the 
second elastic response region and the second yield occurs at a higher force. This is indicative of strain hardening at 
the atomic scale resulting in an increase in the yield modulus and strength. 



1. First Yield : Defect Nucleation 



To gain insight into this behavior, the evolution of the stress profiles in the contact region during indentation are 
analyzed. Figures || and [| show the von Mises shear stress ^fJ-i and hydrostatic pressure p profiles in the region 
directly beneath the indentcr at various stages of indentation marked (a) — (d) in Fig. ^. Figures |(a), (b) and |(a), 
(b) show that as the indentation proceeds from (a) to (&), an increase in ^/J^, a measure of the elastic stored energy, 
substantiates the elastic response seen in the force curve in this regime. At point (b) the elastic stress reach a threshold 
beyond which plastic deformation occurs that partially relieves and dissipates the elastic energy from the surface as 
seen in ||(c) . This behavior is in accordance with the von Mises criterion,t2l which suggests a critical value for VJ2 for 
the onset of plastic activity. Upon further indentation from (c) to (d) \[7i increases again implying another elastic 
response regime, which culminates in the second yield event. 

To study the nature of plastic deformation and characterize the defect structures nucleated, the deformed regions 
are identified by the slip vector soi , where represents the initial undeformed state and 1 represents the state after 
the first yield event. Three snapshots of the deformed region at various stages of defect nucleation between (b) and (c) 
are shown in Fig. || to illustrate the evolution of the defect structure. From the slipped atoms shown in Fig. || it is seen 
that dislocation loops nucleate on the four {111} planes at the surface and extend into the solid. These dislocation 
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FIG. 2: Force versus distance curve during initial stages 
of indentation, (a) — (b) elastic response; (b) onset of first 
yield; (b) — (c) first plastic yield event; (c) — (d) second elastic 
response at a higher force and with a higher slope indicating 
strain hardening. 1 and 2 represent the first and second yield 
events. 




FIG. 3: Contour plots of the atomic von Mises shear stress \f~J2 in the indented region at four stages of indentation marked 
(a) — (d) in Fig. |^. The contours are on [001] surface (upper row) just beneath the indenter and [100] surface (lower row) at 
x = 0. O are the slipped atoms that comprise the defect nucleated during the first yield event and • are the undeformed 
atoms. Stress is concentrated at the corners of the contact region. Increase in \J~~J2 from (a) to (b) and (c) to (d) signifies 
elastic responses and a drop from (b) to (c) indicates plastic yield. 



loops grow in size and intersect with the loops on the adjacent planes forming a pyramidal defect structure as seen in 
Fig. ||(c). Figure ^ shows the corresponding slip vectors of the atoms on one of the slip planes, (111). From Fig. ^|(c) 
the magnitude of the slip vector of the atoms on the (111) plane is close to 1.66 A along [112], which is consistent 
with the (112) partial dislocations on {111} planes in gold. The dislocation nucleated on the (111) plane is therefore 
the i [112] Shockley partial, and similarly |[112], ^ [112] and |[112], partials are nucleated on (111), (111) and (111) 
planes, respectively. Thus the defect consists of intersecting intrinsic stacking faults on the four {111} planes, which 
intersect the (001) surface with four fold symmetry This pyramidal defect structure is in excellent agreement with 
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FIG. 4: Contour plots of the atomic hydrostatic pressure p in the indented region at four stages of indentation marked 
(a) — (d) in Fig. |^. The contours are on [001] surface (upper row) just beneath the indenter and [100] surface (lower row) at 
x = 0. O are the slipped atoms that comprise the defect nucleated during the first yield event and • are the undeformed 
atoms. Stress is concentrated at the corners of the contact region. After the first yield (c), a compressive strain of 0.052 in 
the defect gives rise to increased pressure of the order of 10 — 15 GPa. 
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FIG. 5: Snapshots of deformed region depicting the evolution of the dislocation structures nucleated during the first yield 
event ((b) - (c) in Fig. |). Greyscale represents |sqi|. 

the experimentally observed pcrniajicnt deformation structures^ 

From continuum elastic theory^ of indentation by a rigid flat frictionless punch, similar to the atomic indenter 
used in this study, the stresses reach a theoretically infinite value at the edges of the indenter. This observation, at 
the atomic scale, is validated by the large concentrations piliitrcsses at the periphery of the contact region as seen 
in Figs. ||(b) and |](b), in agreement with previous studies .HtlJ Surface nuclcation of partial dislocations in areas of 
enhanced strepjiS||*asil supported experimentallyETJEJ and theoretical models of surface dislocation nucleation at stress 
concentratorsE30E-iE3 are well established. The corners of the contact region act as stress concentrators and serve as 
the sources of nucleation of Shockley partials on the surface. Under the influence of the stresses, these partials glide on 
the {111} planes, the dominant slip planes in gold, forming intrinsic stacking faults as seen in Fig. ^. This slip results in 
the flow of part of the elastic energy from the contact surface to the sheared surfaces and is seen as an increase in \J~J2 
in the region of slip in Fig. |^(c). The strain undergone by the atoms in the deformed region during defect nucleation is 
quantified by the mean Lagrangian strain tensor in Voigt notation, E = [0.068 0.064 —0.166 0.007 —0.002 0.02]. 
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FIG. 6: Slip vector (soi) maps on (111) plane corresponding to the snapshots in Fig. |5|. O represent the atoms of the slipped 
plane and □ represent the atoms of the unslipped plane adjacent to the slipped region. 



Thus the observed strain is in (100) directions and gives a volumetric strain of -0.052 in the pyramidal defect, giving 
rise to higher pressures of the order of 10 — 15 GPa (with a bulk modulus of 207 GPa) above the ambient pressures 
as seen in Fig. |](c). 

In order to understand why this particular slip system {111}(112) has been activated, we study the resolved shear 
stresses (RSS), as the direction of slip has been associated with maximum RSS. Fig. || shows the distribution of 
T (in)[ioi] i T (in)[oii] 3 pfi r (in)[ii2] on t ne (111) slip plane along [101], [011] and [112], respectively. These are the 
favored slip directionscJ on (111) plane for the current stress state (<r « o zz ). It is seen that even though t^^^oi] and 
T (m)[oiI] have higher concentrations near the surface than T/m^nm, the slip occurs alone [112]. This incongruency 
of slip occurring in a lower RSS direction is consistent with the findings of other. researchersO and can be explained by 
the concept of generalized stacking fault energy (GSF) 7, introduced by Vitek.EZrn3 GSF is the energy per unit area 
of a fault plane created by the rigid slip of one half of a perfect lattice relative to the other along a slip plane in a 
general slip direction. Fig. shows the unrelaxed GSF in the dominant slip directions [110] asd [112]. 

The maximum value of 7 in the direction of slip, called the unstable stacking energy 7 us ,cJ is the energy barrier 
to be overcome during slip. It is seen from Fig. |7] that j us displays a strong directional dependence and slip along 
[112] has a lower energy barrier than [110] and is thus more favorable as observed in this study. However, 7 is a 
static quantity and is not an appropriate measure to describe the dynamics of slip. A better quantity would be the 
theoretical shear stress required to initiate and maintain the slip along the slip direction. This shear stress is given by 



th _ djn (f ,, 

T {«}(r> - -jfr W 

where, 7„ is the GSF on a slip plane with normal n and r is the displacement vector. A plot of this stress for slip 
on {111} plane along (110) and (112) directions is shown in Fig. |7| The maximum value of r th is the resolved shear 
stress that is required to complete the, slip along the particular direction. This is the critical resolved shear stress 
(CRSS) t c referred to in Schmid lawEJ as stated above. From Figs. and || Tfiiiinjm i s 2-34 GPa and is smaller 
than the observed T(m)[ii2] values. On the other hand T^m^no) has a value of 8.88 GPa and is much higher than 
the observed Tm;n rugi values. Thus, the observed slip direction is [112] rather than [011] or [101]. Second and higher 
derivatives of 7 could be used to further refine the dynamics of slip, but it is beyond the scope of this paper and for 
the present study r would suffice. 

The smallest of the directionally dependent r c represents the ideal shear strength of the crystal is 2.34 GPa in 
excellent agreement with the experimental estimates of 1.5 to 2.0 GPa.oO The high r values of 5 GPa, in Fig. ||, 
greater than the theoretical estimate of the ideal shear strength might seem out of order, but it needs to be clarified 
that the theoretical value is based on a block like shear and such instantaneous rigid slip cannot be expected during 
the actual nucleation and propagation of slip. These r values are also quite high compared to the experimental CRSS 
because the values in Fig. |8J are highly localized and are inaccessible to experimental investigations. The shear stress 
values deduced from experiments represent the mean value of the shear stresses in the region local to indentation and 
a plot of such a mean of r is shown in Fig. ^[ It can be seen that these values reach a maximum and drop abruptly 
at the first yield point validating the manifestation of Schmid law at the atomic scale. The maximum RSS on (111) 
is along [112] and reaches a value of 1.95 GPa and on the other {111} planes the maximum RSS are in the range of 
1.8 to 2.3 GPa. These values agree exceptionally well with the experimental estimates of 2 GPa. 

The plastic strains produced byr-indentation are complex and activation of multiple slip systems is necessary to 
accommodate these general yields.E2l Gxoves and Kellytil have predicted the active slip systems by calculating the 
strain produced by a given slip systemEj and identifying the systems that contribute to the observed strain. Such 
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FIG. 7: Generalized stacking fault energies 7 (circles) 
and theoretical shear stresses r" 1 (squares) for slip systems 
{111}(112) (open) and {111}(110> (solid). Lower r ( c 112) indi- 
cates a more energetically favorable {111}(112) slip system. 




FIG. 8: Resolved shear stresses on (111) plane, just before the first yield point, along favored directions of slip: (a) [Oil]; (b) 
[101]; (c) [112]. O are the slipped atoms that constitute the stacking fault and • are the atoms in the undeformed region. In 
some regions r (111)[11 2] is smaller than r (111)[10i] and r (111)[0lI] . 



2.5 




1 0.75 0.5 0.25 o -0.25 -0.5 

Z S ep <A> 

FIG. 9: Mean resolved shear stresses in the deformed region 
for slip systems: □ (111) [011; O (111)[101]; O (111)[112] ; 
* (111) [112]; V (HI) [112]; A (111)[112]. The resolved shear 
stresses reach a threshold at the yield point. The maximum 
resolved shear stress is in the range of 1.8 — 2.3 GPa. 
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(a) (b) (c) 




FIG. 10: Various stages in the formation of the pyramidal defect structure. Shaded regions represent intrinsic stacking faults, 
(a) nucleation and glide of partials; (b) dislocation loop growth and zipping resulting in sessile stair-rods; (c) pyramidal defect 
structure. 

a geometrical analysisEl for compression in [001], which is the observed stress state just before the first yield point, 
predicts activation of slip on the four {111} planes resulting in plastic strains in (100). These predicted slip planes 
and strains are identical to those observed at the atomic scale in this study. 

In the event of the presence of multiple scts.of. independent slip systems capable of producing the required strain, as 
is the case with f.c.c. crystals, Bishop and HillacJ proposed a stress criterion for yielding that requires the attainment 
of CRSS on the active slip systems, without exceeding CRSS on the inactive systems. The observed slip directions 
conform to the above criterion with the shear stresses reaching their critical values in (112), but not in (110). 



2. Defect Nucleation Mechanism 

Based on the results and discussion presented above, we propose a three step mechanism for the formation of the 
pyramidal defect during indentation of Au (001). It is convenient for the following discussion to use Thompson's 
notation for Burgers vectors and planes and refer to Fig. |l^ for an illustration of the mechanism. 

Dislocation nucleation: Surface indentation of Au (001) with an indenter results in large concentration of stresses at 
the corners of the contact region. These stress concentrators, where the RSS on the {111} planes reach the CRSS, act 
as the sources for surface nucleation of Shockley partials (112)) Ba, A/3, 7D and SC on the slip planes (a), (b), (c) 
and (d) respectively. 

Dislocation glide: These partials, under the influence of the external stress due to indentation, glide away from the 
surface forming intrinsic faults on the slip planes. 

Dislocation reaction: As the dislocation loops grow, the partials attract each other in pairs and zip to form sessile 
stair-rods along AC, AD, BD and BC according to the following reactions: 



SC + (3A = S(3/CA 
A/3 + D7 = AD//3 7 
7D + aB = 7Q'/DB 

Ba + C5 = BC/aS (7) 



In vector notation the energetically favorable reactions arc of the type 



i[112]-|[112]=|[100] (8) 

The final defect, therefore, consists of a pyramid of intrinsic stacking faults on {111} planes, which intersect the 
(001) surface with a four fold symmetry, and the (011) edges of the pyramid consist of low energy sessile stair-rod 
dislocations. These sessile stair-rods act as barriers to further glide giving rise to the observed strain hardening during 
indentation beyond the first yield. 
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FIG. 11: Deformed region illustrating the dislocation locks 
formed during the second yield. O represent the atoms that 
slipped during the first yield (|soi| > 0) and O are the atoms 
that underwent slip during the second yield event (|si2 1 > 0). 
Dislocation loops extend beyond the adjacent faces of defect 
structure forming dislocation locks. 

To examine the dependence of the defect structures on the orientation of the indenter with respect to the crystal- 
lographic axes of the gold substrate, the calculations were repeated with the edges of the indenting face in (100) and 
(010) directions. The defect produced was similar to the pyramidal structure seen above where the high stresses at 
the corners of contact region nucleate partials on the four {111} planes. From the physics of the defect nucleation 
presented above, the above mechanism can be generalized to indentation of other {111} and {110} surfaces as well. It 
can be deduced that the indentation of 011} surface will produce a tetrahedral defect structure displaying the three- 
fold symmetry observed in experiments. □ And similarly the hexagonal nature of the indent produced by indenting the 
{110} surface is consistent with the proposed mechanism. Thus, the geometry of the defect structures is independent 
of the indenter orientation, but is a characteristic of the crystallography of the surface of the indented crystal. 

3. Second Yield : Dislocation Locks 

Further indentation results in a second yield, at which point dislocation loops are nucleated on the slip planes 
outside the defect as shown in Fig. O. Figure [l2| shows the slip vectors s 12 , of the atoms on the (111) plane dislocated 
during the second yield and the dislocation loops extending beyond the stair-rods can be seen. The contour plots 
of RSS on (111)(112) along with the slip vectors are shown in Fig. [l3]. It is seen that the activated slip direction 
is not along the maximum RSS direction, but along the direction in which the RSS has reached the critical value. 
This observation further corroborates the discussion presented above. Thus, deformation results through a sequence 
of elastic and plastic responses, with the elastic responses culminating in plastic events. 

B. Retraction 

When the indenter is retracted after the first yield point, the force curve retraces the indentation path at small 
displacements. This suggests that the defect nucleated at the first yield point has disappeared and the substrate has 
recovered is original undeformed state upon retraction. However, the force curve during retraction after the second 
yield point signifies permanent deformation. 

Upon retracting the indenter, the external stress vanishes and the internal compressive stress dominates the dy- 
namics of deformation. It is seen above that the interior of the pyramidal defect is under enormous compressive 
stress due to the strain imposed during the first yield. This compressive stress when resolved onto the slip planes is 
opposite in direction to the RSS during indentation. A restoring force is induced that is strong enough to effect the 
unzipping of the stair-rods into their constituent partials, which then glide toward the surface healing the stacking 
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FIG. 12: Slip vector S12 of the atoms on (111) plane. O rep- 
resent the atoms of the slipped plane and □ represent the 
atoms of the unslipped plane adjacent to the slipped region. 
The slip is along the direction in which the resolved shear 
stress has attained the critical value. 
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FIG. 13: Resolved shear stresses on (111) plane, at the second yield point, along the direction: (a) [121] ; (b) [211] ; (c) [112]. 
O represent the atoms with slip vectors coinciding with the RSS direction; O represent other slipped atoms and • represent 
the atoms in the undeformed region. The observed slip is not necessarily in direction of the maximum resolved shear stress. 



fault along the way. Eventually, the defect disappears and the substrate recovers its original configuration with no 
residual deformation as seen in Fig. ^|. 

However, after the second yield, the dislocation loops extend beyond the stair-rods leading to dislocation locking. 
And upon retracting the indenter after the second yield, the aforementioned restoring forces are not strong enough 
to unlock the locked structure, thus giving rise to permanent plastic deformation observed in Fig. @. 



IV. CONCLUDING REMARKS 



We have investigated the atomistic mechanisms of plastic deformation during nanoindentation of Au (001) surface 
with a noninteracting indenter. A recently developed slip vector analysis has been employed to identify the defect 
structures formed during initial plastic yield. During indentation, the accumulated elastic energy in the indented 
region is partially relieved by the nuclcation of a pyramidal defect structure. The defect is formed by the surface 
nucleation of Shocklcy partials on the four {111} slip planes at the periphery of the contact region. These partials 
glide away from the surface creating stacking faults that grow in size and intersect with those on the adjacent planes. 
At the intersections, the partials zip to form sessile stair-rods which contribute to the strain hardening observed after 
the first yield. The observed slip is in the most energetically favorable direction, which corresponds to the direction 
in which the RSS has reached the critical value and is not necessarily the maximum value. The CRSS estimated in 
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this study is in the range of 1.8 — 2.3 GPa in excellent agreement with the experimental estimates. 

Upon retracting the indenter after the first yield, the pressure due to the compressive strain in the defect induces 
restoring forces that heal the plastic deformation. Further indentation results in a second yield that causes the 
dislocation loops to extend beyond the stair-rods forming dislocation locks. The unlocking forces of these structures 
is greater than the internal restoring forces active during indenter retraction and thus effect permanent deformation 
after the second yield. 

We proposed a three step mechanism based on dislocation theory that elucidates the physics behind the formation 
of the observed defect structures during gold nanoindentation. According to this mechanism, the defects produced 
depend on the crystallography of the indented surface as seen in experiments. 
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